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ABSTRACT 


Time  harmonic  modal  electromagnetic  fields  in  two-medium  half  spaces 
are  investigated.  For  practical  and  numerical  considerations,  the  primary 
sources  of  the  modal  fields  are  chosen  to  be  the  spherical  multipoles,  and 
the  potential  vectors  are  z-directed.  It  is  shown  that  modal  fields  of 
such  combination  are  not  able  to  represent  a  conventional  spherical  modal 
field.  The  horizontal  rotating  potentials  are  added  to  ensure  proper 
representation  and  fast  convergence. 

The  recurrence  relations  which  transform  the  spherical  Hankel -Legendre 
functions  into  the  Fourier-Bessel  integrals  are  derived.  The  secondary 
fields  of  the  Sommerfeld's  type  are  obtained  for  all  spherical  multipole 
sources,  and  the  added  horizontally  rotating  potentials.  The  combination 
of  the  modal  fields  are  capable  of  representing  arbitrary  electromagnetic 
fields  resulting  from  radiation  and  scattering  problems. 

Results  are  given  for  scattered  surface  fields  of  a  buried  finite 


dielectric  cylinder. 


I .  INTRODUCTION 

Modal  expansions  of  electromagnetic  fields  are  fre¬ 
quently  a  mathematical  necessity  in  classical  solutions 
of  electromagnetic  problems.  Recent  development  of  the 
unimoment  method  of  computation  has  also  made  extensive 
use  of  modal  fields,  which  satisfy  the  radiation  con¬ 
ditions  outside  spherical  surfaces  [1] ,  [2] .  Indeed, 
the  application  of  the  modal  fields  is  essential  in 
terminating  the  finite  difference  or  the  finite  element 
equations  at  a  spherical  surface.  It  is  quite  evident 
that  the  unimoment  method  can  readily  be  applied  to 
scattering  or  radiation  problems  involving  air-ground 
half-spaces,  if  the  modal  fields  satisfying  the  air-ground 
interface  boundary  conditions  are  available. 

The  integrals  which  represent  the  potentials  of  a 
dipole  near  an  air-ground  interface  are  known  as  Sommerfeld 
integrals,  which  were  first  derived  by  Sommerfeld  in  1909 
{3]  and  1926  [4] .  They  have  been  the  foundations  of 
research  in  electromagnetic  radiation  by  sources  near  a 
lossy  ground  ever  since.  The  objective  of  this  paper  is 
to  investigate  the  generalization  of  Sommerfeld  integrals 
for  the  modal  fields  so  that  they  can  represent  arbitrary 
electromagnetic  fields  near  an  air-ground  interface. 
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The  generalization  of  Sommerfeld  integrals  has  been 


given  by  C.T.  Tai  [5] ,  whose  systematic  approach  to 
Dyadic  Green's  function  very  elegantly  arrives  at  modal 
fields  of  Sommerfeld  type.  However,  Tai's  generalization 
entails  cylindrical  sources,  i.e.,  each  of  the  exterior 
modes  is  singular  along  the  z-axis.  Using  Tai's  modal 
fields  to  terminate  the  finite  element  mesh,  one  would 
have  to  use  an  infinite  cylinder  as  a  terminating  surface, 
such  as  shown  in  Fig.  la.  The  matching  points  between 
the  finite  element  solutions  and  the  modal  fields  would 
spread  fairly  far  along  the  cylinder,  and  the  numerical 
objective  would  be  finding  a  continuous  function  f(\) 
for  the  values  of  X  along  the  integration  path.  The 
above  particulars  are  not  attractive  features  for  compu¬ 
tational  purposes.  The  ideal  matching  surface  should  be 
finite  and  the  ideal  numerical  objective  should  be  a 
discrete  set  of  coefficients.  In  this  paper  we  shall 
attempt  to  generalize  Sommerfeld  integrals  based  on 
spherical  multipole  expansions.  The  numerical  application 
of  such  generalized  Sommerfeld  integrals  should  result  in 
spherical  matching  surfaces,  as  shown  in  Fig.  lb,  and 
discrete  numerical  objectives. 

The  fields  concerned  in  the  following  discussions 
are  all  assumed  to  have  the  exp  (jcot)  harmonic  time 
dependence . 


_ 


II.  MULTIPOLE  EXPANSIONS  AND  VERTICAL  VECTOR  POTENTIALS 


The  first  step  in  deriving  Sommerfeld  integrals  is 
to  decide  what  type  of  primary  sources  are  to  be  used. 
The  electromagnetic  fields  in  a  homogeneous  isotropic 
medium  can  usually  be  derived  from  the  electric  and 
magnetic  vector  potentials,  Ag  and  Am  as  follows 


E  =  V  x  A  +  —  7  x  7  x  A 

m  30)£  e 


^  1  -*• 
H  =  7  x  A  -  “t— — *  7  x  7  x  A_ 


(la) 

(lb) 


e  3ojy  m 

where  e  and  u  are  the  permittivity  and  the  permeability 
of  the  medium,  w  is  the  angular  frequency  and  the  potentials 


A  and  A  both  satisfy  the  vector  Helmholtz  equations.  The  above 
e  m 


potential  representations  are  complete  in  homogenous  source  free 
medium  where  the  fields  derived  from  the  gradient  of  a  potential 
vanishes  [6]. 

Our  intention  to  obtain  spherical  matching  surfaces 
in  numerical  application  limits  our  option  to  spherical 
multipoles,  i.e.,  the  primary  fields  will  be  derived  from 


=  z 


E  E 


00  /a. 


m,n 


i*h^  (kr)  Pm(cos  0)  e^m<^ 
k  n  n 


m, 


m=-«  n=|m|\B 


m,  ny 


(2) 


where  k  =  w/pT  is  the  wave  number  and  z  is  a  unit  vector 
normal  to  the  ground  plane.  It  is  noticed  that  the  vector 
potentials  are  z-directed  rather  than  the  conventional 
radial  vector.  The  direction  of  the  potential  vector  is 
dictated  by  the  need  to  satisfy  the  continuity  conditions 

of  the  planar  air  ground  interface  for  the  total  fields. 
The  use  of  the  conventional  spherical  potentials  using 
radial  vectors  does  not  lead  to  expressions  which  are 
compatable with  air-earth  interface  field  matching. 

This  particular  combination  of  spherical  harmonics  and 


z-directed  potential  vectors  is  not  new.  The  complete 
expression  has  been  given  by  Tai  [5  Appendix  B] ;  however,  the 
modal  fields  derived  from  it  have  never  been  applied  to  solve 
specific  problems.  For  the  convenience  of  the  further 

discussions,  the  fields  derived  from(2)  will  be  named  as 

(z)  (z)  -*■  ■* 

TM  and  TE  modes  for  A  and  A  respectively. 
m,n  m,n  e  mr  1 

The  modal  electromagnetic  fields  derived  from  (2)  are 

denoted  by  g™(2)  ^TM(z)  an(^  ^TE(z)  ^TE(z)  >phe  definitions 

J  m,n  m,n  m.n  m,n 

of  these  notations  are  listed  in  Table  1. 

One  of  the  most  important  questions  about  the  modal 
fields  derived  from  (2)  is  whether  the  expansions  can  rep¬ 
resent  arbitrary  radiating  fields.  The  follow¬ 
ing  derivations  will  show  that  the  modal  expansions  of  a 

class  of  radiating  fields  by  TM  and  TE'  _  result  m 

^  J  m,n  m,n 

diverging  series.  Such  a  class  of  radiating  fields  include 
the  conventional  spherical  vector  wave  modes.  It  is  found 
that  proper  additional  terms  may  be  added  to  (2)  to  make 
the  modal  expansion  series  convergent. 

For  clarity  in  the  ensuing  discussions,  we  should  examine 
the  conventional  curvilinear  spherical  vector  waves  obtained 
from 


co  co  /  a 

II 

m=-°°  n=)m|  \b 


hn2)  (kr)  pn (cos  6)  ejm<^ 


where  r  is  a  radial  vector.  The 


se  expansions  were  proved 


— •v- 
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to  be  complete  in  the  homogenous  source  free  region  where 
the  divergence  of  E  and  H  vanishes  [6]  and  have  been  applied 
to  numerous  practical  applications.  We  shall  name  the  ex- 

/  „  \  /  „  V  ^  ^ 

pansion  modes  of  (3)  as  TH  '  and  TE '  for  A  and  A 
r  m,n  m,n  e  m 

respectively.  The  modal  fields  are  denoted  by  e^M^r^ , 

£TM(r)  and  |TE(r)  £TE(r)as  shown  in  Table  1. 

m,n  m,n  m,n 

We  first  study  the  singularities  of  TE^,  TM^r^ 

3  m,n  m,n 

and  TE  ^ z ^  ,  TM^  fields  as  r  +  0.  The  spherical  Hankel 
m,n  m,n 

function  and  its  first  derivative  have  singularities  at 
r  -*■  0 

lim  h^2) (kr)  -  (kr) ~n 
r-0 


lim 

r+0 


d  h^2) (kr) 
d  kr 


-j (kr)~(n+1) 


Hence  the  modal  electric  fields  have  the  singularities 
as  follows: 


lim 

r-*0 


lim 

r-0 


-TE (z) 
em,n 


-TM(z) 

em,n 


r- (n+1) 
r-(n+2) 


lim 

r-0 


±TE(r) . _ 
m,n 


r 


-n 


lim 

r-0 


gTM(r) 
em,  n 


r~ (n+1) 
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The  expansions  of  different  numbers  of  m  are 

decoupled  because  the  Fourier  series  modes  are  mutually 

orthogonal.  For  a  particular  value  of  m,  we  found  that 
( z )  ( z ) 

as  r-*-0  TE  and  TM  have  a  singularity  sequence  of 
m,n  m,n 

r"(|m|+D,  r-(lml+2),  ...  whereas  TE^  and  TM^  have 

r  ,  r  ...  Therefore  the  lowest  singularity 

term  r  ^  is  missing  for  the  TE^,  TM^  modes.  This 

m,n  m,n 

suggests  that  additional  terms  of  proper  singularity 

f  2  )  (  Z  ) 

must  be  added  to  TE  and  TM  modes  to  represent 

m, n  m ,  n 

arbitrary  radiating  fields. 

iz)  (z\ 

To  show  that  the  TE  and  TM  modes  do  not  yield 

m,n  m,n  J 

converging  series  when  expanding  a  certain  class  of 
radiating  fields,  let  us  take  a  simple  example  of  the 
fields  radiated  by  a  horizontal  rotating  dipole.  The 
vector  potentials  of  the  dipole  can  be  written  as 


Ae  =  (x  +  jy)  h^2) (kr) 


The  electromagnetic  fields  generated  by  (4)  are  found 

(it) 

to  be  exactly  the  same  as  those  of  the  TM^  ^  mode.  We 

shall  find  the  coefficients  a  ,  8  in  the  following 

n  n 


expansion 


n=lL 


T 


modes  are  ortho- 


It  should  be  noted  that  neither  e™*z'nor 

m ,  n  m,n 

gonal  on  the  surface  of  a  sphere. 

By  using  the  relations  given  by  Tai  [5,  p.  228],  such  as 

-*TE  ( z )  _  •  1  1  ^TM(r)  _ n _  +TE(r)  +  n+1  £TE(r) 

el,n  3  n (n+1) k  el,n  (2n+l) (n+1)  l,n+l  (2n+l)n  l,n-l 


-►TM ( z )  _  .  1  .  +TE(r) ,  _J2 _  ^TM(r)  nfl _  ^TM(r) 

l,n  “  3  n (n+1)  K  el,n  (2n+l) (n+1)  1, n+1  (2n+l)n  l,n-l 

and  substituting  (5)  into  (6a)  and  (6b),  one  can  solve  an  and 

8  via  a  set  of  linear  eauations  involving  a  tri-diagonal  matrix 
n 

The  results  are 


an  = 


(-1) 


( 2n+l) 
n (n+1) 


if  n  is  odd 

if  n  is  even 
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n 


(2n+l) 
n (n+1) 


if  n  is  even 


if  n  is  odd 


Examining  (6a)  and  (6b)  one  finds  that  an  and  8n  have 
asymptotic  behavior  of  0(n  ^)  for  large  orders  of  n.  The 
asymptotic  form  of  the  spherical  Hankel  function  for  large 


(6a) 

(6b) 

(7a) 

(7b) 


n  [7,  eq.  37]  is 


^  (2)  „__4  _  ...  /  1  f  2(n+h) 

lim  hn  (kr)  -  lyicr)  (n+L)  •  (  e. (kr) 


(n+4) 


where  e  is  the  base  number  of  the  natural  logarithm.  Com¬ 
bining  the  slowly  convergence  of  a and  Bn  and  the  rapidly 
increasing  of  h  ^  (kr)  at  (n  +  ?s)n,  it  is 
thus  shown  that  the  expansion  of  (5)  diverges  for  all  kr. 

Similar  conclusions  can  also  be  obtained  for  other  TM_[r^ 

m,  n 

modes. 

The  fact  that  TM.iz^  and  TE^  yield  diverging  series 
m,n  m,n  33 

( r ) 

when  expanding  the  field  of  a  rotating  dipole  and  the  TMm  n 

(  3T  ) 

and  TE  modal  fields  has  made  the  expansion  of  (2)  un- 
m,  n 

usable  for  general  applications  in  scattering  problems.  To 
remedy  this  situation,  additional  terms  of  proper  orders  of 


singularities  must  be  included  in  (2) . 


It  was  found  that 


the  field  of  a  rotating  dipole  given  by  (4)  can  provide  the 
singularity  needed  for  m=l.  Hence,  we  suggest  adding  two 
circularly  polarized  vector  potentials,  e.g., 


=  (x± jy)  £  YifflUh'2)  (kr)?”  <*>«e)e±*m'*  <8> 


A  A  , 

where  x  and  y  are  unit  horizontal  vectors  on  the  rectangular 
coordinates.  These  rotating  modes  will  be  named  R™+ (m+i) 
and  RTE±^m+1j  for  Afi  and  Am  respectively.  The  modal 
fields  will  be  denoted  by  ej™+1)  ,  Sj™+1)  and  ' 
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"*RTE 

h+(m+i).  Indeed,  these  rotating  modal  fields  are  similar 

to  the  lowest  order  modes  for  each  m  of  the  conventional 

spherical  vector  waves.  That  is,  they  are  similar  to  the 

modal  fields  of  and  TE_fr^ 

±m,m  ±m,m 

The  proper  general  solution  of  Maxwell's  equations  out¬ 
side  a  sphere  in  terms  of  (2)  and  (8)  is  thus 


oo  oo 


J-  E  l 


m=_oo  n=  m 


*  TE(z) 
"m,n 

m,n  '  £TE(Z) 
hm,n 


m=l 


*E 

m=l 
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-m 


*RTE 

'm 


"RTE 
1m 


+  n 
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►RTE 

J-m 


►RTE 

'-m 


+  n 


+  B. 


-m 


nvn 


TTM(z) 
m,  n 


►TM(z) 
‘ra,n  4  -* 


-m 


where 


a  .8  ,  y  ,  and  n_ 

m,n  '  pm,  n  'm  ™ 


(9) 


are  the  expansion  coefficients.  The  modal  fields  of  the 

expansion  are  defined  in  Table  1. 

It  can  be  shown  [8]  that  any  one  of  the  conventional 

spherical  vector  modal  fields  can  be  expanded  in  a  finite 

number  of  terms  by  using  (9) .  The  convergence  rate  of  (9) 

is  the  same  as  that  using  TM(r*  and  TE ^  modes  when 

m,n  m,n 

solving  the  same  problem. 
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Table  1.  The  Definition  of  the  Primary  Modal  Fields 


MODES 

VECTOR  POTENTIALS 

MODAL  FIELDS 

te|z) 

tm,  n 

A  *  2  r  h*2^  (kr)P™(cos6)e±^m<*' 
m  k  n  n 

3“(*>-  PxA 
±m,n  m 

RTE(z)  ■  1_  VxVxA 

n±m,n  3  kZ  VXVXAm 

mm 

Ag  =  z  ^  h^2*  (kr)  P™  (cos9)  e^m<^ 

|TM(z)=  ,Z  Vx7x^  =  _z2jjTE<*> 
eim,n  ->k  VXVXAe  ±m,n 

K^M(Z)=  VxA  =  £JE(z) 

±m,n  e  ±m,n 

nTEim 

Sm  *  t*-™'  E  hm-l  (!tr) 

P™-f  <cose)e!3(”-1,» 
m-1 

-  7»km 

gRTE  •  1 VxVxA 

n±m  -1  kZ  VXV  m 

RTM±m 

■  (Xijy)  I  ^5.  (kr) 

Ci  (cos9)e  C”"11  * 

£R™  =  -ji  VxVxA  =  -Z2hRTE 
e±m  ^k  xve  in 

hR™  =  VxA  =  eRTE 
n±m  VXAe  e±m 

T>f  ^  ^  ^ 

'~±m, n 

A  =  r  h(2)  (kr)  Pm (cos6)  e*^™^ 
m  n  n 

7** 

±m  ,  n  m 

jTE (r)  ,1_  v  7  ^ 

ft±m,n  DkZ  VXVXAm 

TM^r) 

tmr 

Ae  =  r  hj2)  (kr)p"(cos0)e±jm<l> 

S™'?-  jl  7*7*A  -  -22Slf<7> 
±m,n  Jk  e  ±mfn 

jjTM(r)  7  -  =  -TE(r) 

±m,n  XAe  ±m,n 

NOTE:  Z  =  /  u/e 
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III.  THE  FOURIER- BESSEL  INTEGRALS  OF  h(2) (kr)  Pm  (cos6) 

- — ■ — — - - - n - n - 

The  multipole  expansions  discussed  in  the  last 
section  are  for  the  primary  fields,  which  are  valid  in  an 
infinite  uniform  space.  To  find  the  fields  of  the  multi¬ 
poles  near  an  air-ground  interface,  the  multipole  fields 
must  be  expressed  in  cylindrical  harmonics  so  that  the 
.secondary  waves  can  be  derived.  We  shall  derive  the 
Fourier-Bessel  integrals  for  the  spherical  Hankel- 
Legendre  functions  in  this  section. 

Since  both  h^(kr)  Pm  (cos  0)  e_^m<^  and 
n  n 

— 2 - T 

*  k  ^  e^1"^  are  solutions  of  the  same  scalar 
m ' 

Helmholtz  equation 

7  2<l>  +  k2t|/  =  0  (10) 


one  may  represent  the  spherical  harmonics  by  the  superposition 
of  a  complete  cylindrical  eigenfunction  [6,9].  That  is 

.  CO 

h(2)  (kr)  P*  (cos  9)  =  [sgn  ( z)  ] n_m  J  n  (A)  Jm(AP)  e'ulzldA  (11) 
n  ...  0 

I  2  2 

where  u  ="\J A  -  k  ,  Real(u)  *  0,  and  sgn  (z)  =  1  for  positive 
z  and  -1  for  negative  z.  The  function  fm  n(A)  is  an  amplitude 
function  for  the  transformation.  The  first  special  case  of  the 

Fourier-Bessel  transform  in  (11)  is  for  m=n=0  of  a  dipole, .as  shown 
by  Sommerfeld  in  [9]. 

h02>  ,krl  *  1  TT  =  /  3  >T5  VXo>  e"u|zl  ax  <12) 

0 
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Hence  the  function  ^(X)  is 

‘o,o(X>  -  5  m  (13) 

The  functions  for  other  orders  of  m  and  n  can  be  obtained 
by  recurrence  relations. 

The  following  two  recurrence  relations  are  proved  in 
Appendices  A  and  3,  respectively. 


Recurrence  Relation  1. 


hm+l  lkt>  Cl  lcOS  91  *  -<2"+1)  Wi  thm2’  (kr)  Pm  (cos  9)1 

+  [h(2)(kr)  P®  (cos  6)] 

Kir  sin  c  m  in 


Recurrence  Relation  2 . 


fhn2)(kr)  pn  (cos  8)] 


=  (n+m) _  h^2^  (kr)  Pm  (cos  0) 

(2n+l)  n-1  n-1  lCOS  ’ 

_  ( n-m+ 1 )  ^(2)  (kr)  Pm  (cos  9) 
( 2n+l)  n+1  (  '  n+1  (  J 


The  above  recurrence  formulas  are  unique  in 
that  they  relate  spherical  modal  potentials  rather  than 
single  variable  special  functions.  The  physical  interpre¬ 
tations  of  these  formulas  are  also  of  interest,  but  we  shall 
not  attempt  to  discuss  these  here. 


It  is  noted  that  the  first  recurrence  relation  raises 
the  order  of  m  ,  and  the  second  raises  n  for  a  fixed 
value  of  m  . 

The  formulation  of  f  , .  , . (A)  is  related  to  f  m  (X) 

m+l ,m+l  m,m  ' 

by  substituting  the  integral  of  (11) into  the  first  recur¬ 
rence  relation  (14).  That  is 


00 

/ 


fm+l,m+l(X)  Jm+l(Xp)  ‘  dX 


=  - (2m+l) 


00 

ft 

J  m, 
0 
00 

/ 


m 


(A) 


3  J__  ( A  p ) 

m  in  w  t  y  \ 

- 51 - -  i —  J  (  Ap) 

3kp  kp  m 


=  ( 2m+ 1 )  /  f  _(A)*  A  J_li  (Ap)  e"U'2'  dA 

*  m m  k  m+ 1 


e-^l  dX 


(16) 


This  relation  holds  for  all  values  of  p  and  z,  hence 


^m+1 ,m+l ^ X^ 


(2m+l) 


A 

k 


f 


(A) 


(17) 


Using  the  initial  formula  of  f^  g(A) 
(171 ,  one  has  a  general  expression  of 


f  mU) 

m,  m 


Pm(0) 

m 


in  (13)  and  the  relation 

f  _(A)  as  follows: 
m  ,m 

(18) 


for  m=0,  1,  2 .  Because  Pm  .  (cos  9)  =  0,  we  have 

m- 1 


,  (A)  5  0 
m,  m-1 


(19) 


Similarly,  by  substituting  the  integral  of  (11), the 
second  recurrence  relation  (15)  leads  to  the  following 


relation: 


Recurrence  Relation  3. 


fm,n+l<X> 


(2n+l)  u 
(n-m+1)  k 


f 

m,  n 


(1)  + 


(n+m) 

(2n+l) 


f 

m,n- 


1 


(20) 


A  complete  formulation  of  fm  n(X)  can  readily  be  obtained  by 
(18)/  (19)  and  the  third  recurrence  relation  (20).  The 
recurrence  relations  ( 17 )  and  (20)  resemble  those  for  the 
associated  Legendre  functions.  We  found  that  f^  n(X)  can 
be  written  in  a  closed  form  to  be 


f 

m,n 


(X) 


.(n-m+1)  X  vm  .  .  u. 
3  ku  n  l~3  k 


(21) 


Although  (21)  may  also  be  obtained  by  other  approaches, 
they  have  not  been  explicitly  shown  in  open  literature. 

The  approach  we  have  presented  here  is  more  direct  and 
computationally  appealing  than  integration  of  plane  waves  in 
the  complex  plane. 

It  is  worth  noting  that  the  Fourier-Bessel  integral 
is  equivalent  to  the  Fourier-Hankel  integral  [9], 


h(2) (kr)  Pm  (cos  9) 
n  n 


(  sgn  (z)  ] 


n-m 


(X)  H^2)  (Xp)  e'u,z 
2  J  m, n  m 


dX 

(22) 


The  integration  paths  and  branch  cuts  in  the  complex  X-plane 
for (ID  and  (22)  are  shown  in  Figure  2.  The  path  is  a 

permissible  integration  path  of  the  Fourier-Bessel  integral 
(21),  and  path  P2  +  Pj_  is  for  the  Fourier-Hankel  integral  (22). 
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IV.  THE  GENERALIZED  SOMMERFELD  INTEGRALS 

Let  the  multipoles  be  at  a  depth  d  underground. 

The  coordinates  in  the  meridional  plane  are  shown  in  Figure  3. 
The  relative  dielectric  constants  of  regions  I  and  II  are 
and  respectively.  We  assume  the  same  magnetic 

permeability  in  both  regions,  although  it  is  straightforward 
to  consider  the  case  of  different  magnetic  permeabilities. 

-*•  T 

The  vector  potentials  m  regions  I  and  II  are  A 
-►I  ■>‘II  ■*‘11 

(or  A  )  and  A  (or  A  ) .  The  fields  in  region  II  are 
erne 

decomposed  into  primary  and  secondary  waves  as  follows: 


A11 

rfPri 

=  A 

+  ASeC 

(23. a) 

m 

m 

m 

A11 

-*Pri 
=  A 

+  ASeC 

(23. b) 

e 

e 

e 

The  primary  waves  are  the  multipole  fields  in  an  infinite 
space,  which  have  been  discussed  in  previous  sections.  The 
secondary  waves  are  the  fields  due  to  the  air-ground  interface. 

Owing  to  their  basic  differences  in  the  formulation, 
we  shall  discuss  the  derivations  for  TE(z),  TM(z) ,  RTE ,  and 
RTM  fields  separately. 

The  Vertical  Magnetic  Multipoles,  TE(z)+m  n  Fields 

The  primary  field  for  TE(z).  in  region  II  are 

±m,  n  J 


1. 


The  sufficient  conditions  for  the  continuity  of  tangential 


electric  and  magnetic  fields  (Ep,  ,  Hp/  and  )  on  the 
interface  between  regions  I  and  II  are 


and 


3A1 

m 

\.II* 

3z 

z=d 

3z 

(29. a) 


(29. b) 


Substituting  (23. a),  (24),  and  (27)  into  (29. a) 
and  ( 2 9 . b) ,  we  have 


OO 

±jm<J)  f  rf 

J  [  m, n 


-u-d 

.  X )  e  +  g  (  X) 


-  h(X)j 


J  (Xp)  dX  =  0 
m 


J  ru 

r\  L 


-u_d 

f  (X)  e  -  u_g(X)  -  u  h  ( X ) 

2  m, n  2  1 


( 30. a) 


Jm(Xp)dX  =  0 


(  3,0. b) 


The  conditions  of  (30.  a)  and  (  3 0.  b)  must  be  true  for  all  <}> 
and  p  ,  hence 

-u  d 

f  (X)  e  £  +  g(X)  -  h  ( X)  =  0  (  31. a) 

m ,  n 

-u_d 

u,f  (X)  e  1  -  u.g(X)  -  u.h(X)  =  0  (31.b) 

2  m, n  2  1 

The  functions  g ( X )  and  h(X)  are  obtained  by  solving 
(  31. a)  and  (  31.b)  . 
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g(X) 


U2  ~  U1 
u2  +  U1 


f 

m,  n 


(X) 


h  (X) 


2u 

- ~r~ -  f  (X) 

U2  +  U1  ra,n 


(32. a) 

(32. b) 


Replacing  g(X)  and  h(X)  in  (26)  and  (27;  by  the  above 

expressions,  we  obtain  the  solution  of  TE(zK  fields 

±m,n 

for  z  >  d  and  z  <  d  as  follows 


l1  =  2  e±jm* 
m 


f 

J  U2  +  m,  n 


-u_d-u1 (z-d) 

(X)  Jm(Xp)  e  dX 

(33. a) 


l11 

m 


=  z  e 


±jmJh(2)(  pm 

n  2  n 

00 

r  u2  ‘  ui 

J  fB|n(i»  Jm(Ap) 


(cos  9) 


-u  (2d-z)  ] 

e  dX  (33. b) 


Taking  m  =  0  and  n  =  0,  one  finds  that  the  solution 
given  in  (3«3.a)  and  (33. b)  for  TE(z)Q  Q  is  exactly  the  same 
as  that  of  the  Sommerfeld  integrals  for  a  vertical  magnetic 
dipole  buried  underground  at  a  depth  of  d.  Note  that  the 
Hertz  potentials,  Trg  and  irm,  as  used  by  Sommerfeld  and  the  vector 
potentials  used  here  differ  by  constant  multipliers  as 


Ao  =  3^, 


A  = 


m 


m 


(  34) 
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2. 


Fields 


The  Vertical  Electric  Multipoles  ™(z)+m/n 
The  general  expressions  for  fields  of  ™(z)+m>n 
regions  I  and  II  are  identical  to  Eqs.  (23),  (24)  and  (25)  by 


replacing  A_  with  A  . 

The  electric  and  magnetic  fields  corresponding  to  the 

TM(z)  vector  potentials  are 

±m,  n  c 


£  =  -ri¬ 

ll  toe 


*2K  -in,  3A  /  . 

f  —  ‘’ViT*4  (k  Ae  + 


3z 


32A 


3z‘ 


H  =  ±p  ^  A  + 
P  e 


3p 


(35. a) 

(35. b) 


The  sufficient  conditions  that  the  tangential  components  of 
electric  and  magnetic  fields  be  continuous  on  the  interface 
between  region  I  and  II  are 


1_ 

z. 


r  Ti 

3A-1 

r  3  a11*] 

e 

1 

e 

3z 

z=d  Z  2 

3z 

z=d 

■  KL 

Jz=d 


('36* a) 


(36. b) 

Enforcing  the  boundary  conditions  at  the  air-ground  interface, 
we  have 

«  r  , 

u„  -u„d  u„  u. 


s ±  jm4>  J 


-u_d  u9 

~  f  ( 1 )  e  -  g(X)  -  -i  h(X) 

£  2  HI  /  n  E  2  E  ^ 


J  (Xp)  dX  =  0 
m 


(37. a) 


00 

J 


-u_d 

f  (X)  e  +  9(^)  -  h(X) 


J  (Xp)  dX  =  0  (  37. b) 

m 
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* 


it1 

e 


it11 

e 


±^m<J>  f  — -  f  (X)  j  (Xp)  e 

J  £1u2  +  e2u1  m,n  m 


=  z  e 


=  2  e 


-u_d-u, (z-d) 


2  “1 


dA 
(40. a) 


(2) (k_r)  Pm  (cos  0) 
2  n 


CO 

I 


e.u_  -  e  u  -u_ (2d-z) 

-±— £ - f  (X)  J  (Xp)  e  1  dX]  (40. b) 

elU2  +  £2U1  m,n  m  J 
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r 


4* 


The  solutions  given  in  (40-a)  and  (40*b)  for  TM(z)Q  Q 
are  exactly  the  same  as  those  of  the  Sommerfeld  integrals  for 
a  vertical  electric  dipole  buried  underground  at  a  depth  d  , 
except  for  the  differences  between  the  vector  potentials  and 
the  Hertz  potentials,  as  given  by  (34). 

3.  The  Rotating  Magnetic  Multipoles ,  RTE+m  Fields 

The  primary  fields  of  RTE+m  for  m  =  1,2,3,...  in 
region  II  are 


±  jy)  h(2>  (Jc0r)pmr?-  (cos  0) 


00 

'  <Mj*>  Vl11'1 


-U- |z I 
e  z  dX 


For  the  same  reason  pointed  out  by  Sommerfeld  for  a  horizontal 
dipole,  it  is  required  that  the  vector  potentials  of  the 
secondary  fields  have  a  z-component  in  order  to  satisfy  the 
boundary  conditions  of  the  air-ground  interface.  Therefore 
we  assume  the  secondary  fields  in  region  II  and  the  total  field 
in  region  I  as 


t  Sec 
A 

m 


(P  ±  j$)  A I  +  A*  z 
m  mz 


(42. a) 


Am  =  (P  t  j$) 


Ate  +  Afc  2 
m  mz 


(42. b) 


The  sufficient  conditions  for  the  continuity  of  the  tangential 
electric  field  at  z  =  d  are 
.Pri 


*3APri  3Ar"l  f  3Afcl 

L  J2=d  l  Jz=d 


(45. a) 


[a'  I  -  [a1] 

L  m2J  z.d  l  ”ZJ  z=c 


(45. b) 


The  conditions  for  the  continuity  of  the  tangential  magnetic 
field  at  z  =  d  are 


[APrl  +  Ar]  =  e .  [  Afel 

L  m  mJz=d  H  mJz=c 


(45. c) 


r  3Ar 

mz 

3Afc  1 
mz 

p3Afc 

m 

(m-D  tl 

[3z 

3z 

i 

II 

N 

Lap 

P  AmJ 

-[ 


z=d 


~.Pri  azvr 

3Am  +  _m  .  (m-1)  APri  _  (m-1) 
3p  3p  pm  p 


Arl 

mJz=c 


(45. d) 


Substituting  (43. a)  and  U3.c)  into  (45. a)  and  (45. c)  ,  one 
can  solve  for  g(X)  and  h(X)  . 


eiu->  ~  £->ui  -u_d 

g(X)  =  - 12  21  f  (A)  e  2 

e^u2  +  ^2U1 


2£2U2 

U  M  \  _  Z  Z  ~ 

e1u2  +  e2u1  m-1, m-1 


(X)  e 


-u2d 


(46. a) 


(46. b) 


The  functions  g  (X)  and  h  (X)  are  derived  from  (45.b) 

Z  Z 

and  (43 .d) . 


(e.  -  e  )  2u_ X 

gz(X)  “  hz(x)  “  (Uj^  +  u2)  (e1u2  +  e2u1)  fm.-l,m-l 


( X)  e 


-u2d 


(  46.  c) 
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The  complete  vector  potentials  of  RTE^  fields  for  z  >  d 
and  z  <  d  are  readily  obtained  as  follows: 


A1  =  <x  +  i?)  e±j(m~1H  f  -  2-2-? -  f  .  .  (X)  J  ,(Xp) 

'x  +  jY'  e  J  £,u.  +  e..u_  m-l,m-l  m-1 


+  ze 


00  ( 
±jm4>  f  _ 

i  (U1 


-u2d-u1 (z-d) 


(e  -e-)2u-X 

L+u2)  (£2U1  +  £1U2)  Jm(Xp) 

-u_d-u. (z-d) 

^  -L  A  \ 


(47. a) 


a”  -  IS!  j?>OV>  Cl  ,cos  81  S±j(”"1H 


+  (* 


f  el_u2  ~  e2ul  f  m 

“  J  e2ul  +  elU2 


-u- (2d-z) 

J  . (Xp)  e  ,  dX 

m-i 


+  iniA  r  Ul  "  £->)  2u-)A 

t  |  .1  2  2 _  f  (• 

J  (u1 + u2) (e2u1 +  e1u2)  m-1, m-1 v 


— U- (2d-z) 

J  (Xp)  e  ^  dX 

m 


(47. b) 


The  solutions  given  in  (47. a)  and  (  47. b)  for  RTE  +  ^ 
are  the  same  as  the  Sommerfeld  integrals  of  a  buried 
horizontal  magnetic  dipole  at  a  depth  d  rotating  counter¬ 
clockwise  (for  RTE+^)  or  clockwise  (for  RTE_^) . 


4. 


Fields 


The  Rotating  Electric  Multipoles,  RTM+m 

The  general  expressions  for  the  primary  fields  of 

RTM+m  for  m  =  1,2,3,...  in  region  II  are  identical  to  Eqs. 

(41),  (42)  and  (43)  by  replacing  A  with  A  and  A  with  A  . 

m  6  mz  6  z 

The  electric  and  magnetic  fields  for  the  RTM  modes  are 


E  = 


-j~  f  fk2A  +  i-is.  -  (m'AI  a!  0  + 

ue  IL  e  3p2  P  3P  p2  ^ 

±j  fk2A  +  h  .  »is=ii  A1  $ 

L  e  p  3p  p  ej 


32A 


ez 


3p3z 


m  9A 
m  mz 

P  3z 


♦ ra2ft- . 

|_3z3p 

(m-1)  3  Ae 

P  3z  _ 

2  + 

t2%z  *  „r]  0 

( 48. a) 

H  = 

3A 

3A 

+  TT  $  ±j 

*3A 

e 

.^P~  "  ' 

±3  ?  Aez  * 

3A 

ez  * 

"  3  P  $ 

(  48-b) 

The  conditions 

for  the  continuity 

of  the  tangential 

electric 

field 

at  z  =  d 

are 

Krt  + 

SAe] 

K] 

(49. a) 

[_3z 

3zJz=d 

L3z  X- 

=d 

[*"]*-*  *  [*-U 


(49. b) 
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The  conditions  for  the  continuity  of  the  tangential  magnetic 
field  at  z  =  d  are 


[*.'"*  *:L  ■  [*:].« 


(49  -c) 


1 

raAr  1 

ez 

_1_ 

fm  t  *1 

ez 

_1_ 

r~,t  1 

aAe  (m-1)  t 

e2 

3z 

z=d  £1 

3z 

_  «■ 

z=d  £1 

3  P  P  e 

,Pri  .Ar 
3A„  3A 

_§ _  +  _® 

3  p  3  p 


z=d 


(m-1)  APri  _  (m-1)  Ar 
p  e  ~  p  e 


z=d 
(  49.  d) 

Enforcing  the  boundary  conditions  at  the  air-ground  interface, 
we  have 


u-  -  u. 

g(X)  =  -2— r~,4  f 


h  ( X)  = 


u2  +.  ux  m-1, m-1 


2U2  £ 

u2  +  u^  m-1, m-1 


(X)  e 


-u2d 


(X)  e 


-u2d 


(  50-a) 


(S^-b) 


The  functions  g  (X)  and  h  (X)  are  obtained  from  (  49.b) 
z  z 

and  (  49.  c)  . 

(£  -  e  )2u-X  -u  d 

*z(X>  -  VX>  *  (tlU;  %  u2)  £m-l,m-l(X)  e 


(  50.c) 


The  complete  vector  potentials  of  RTM+m  fields  for  z  >  d 
and  z  <  d  are  obtained  as  follows: 
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V.  RESULTS 

Utilizing  the  exterior  modal  fields  generated  by  the  Generalized 
Sommerfeld's  Integrals,  we  are  able  to  compute  the  scattered  surface  fields 
of  a  buried  dielectric  finite  cylinder.  The  scattering  configuration  is 
shown  in  Fig.  4.  Figure  5  shows  the  surface  area  where  the  fields  will  be 
computed.  Figures  6-10  displace  the  3-D  and  contour  plots  of  the  scattered 
fields  for  700  MHz  and  1000  MHz. 


VI.  CONCLUSION 


There  are  several  numerical  ways  to  solve  radiation 
or  scattering  of  electromagnetic  fields  near  a  lossy  half 
space;  notably,  the  method  of  moments  [  11]  ,  the  extended 
boundary  condition  method  [  12] ,  and  the  unimoment  method 
[1] .  All  the  methods  in  one  way  or  another  make  use  of 
the  celebrated  Sommerfeld  integrals.  The  problem  may  also  be 
considered  as  a  limiting  case  of  scattering  by  two  non- 
concentric  dielectric  spheres  with  one  imbedded  in  the  other 
[13] .  But  that  special  approach  entails  such  complex 
addition  theorems  and  integral  expressions  that  only  the 
zeroth  order  azimuthal  mode  has  been  obtained. 

In  this  paper,  we  have  presented  the  qroundwork  for  the 
application  of  the  unimoment  method  to  solve  the  scattering 
by  buried  obstacles,  which  requires  the  generalization  of 
the  Sommerfeld  integrals  for  multipole  sources.  For  success¬ 
ful  numerical  applications,  we  find  it  necessary  to  use 
spherical  harmonics  combined  with  z-directed  potential 
vectors.  The  modal  fields  so  obtained,  however,  appear  to 
be  inadequate  in  representing  a  conventional  spherical  mode. 
While  the  situation  is  duly  corrected  by  addition  of  two 
horizontally  rotating  potential  vectors  for  the  lowest  order 
azimuthal  modes,  further  theoretical  investigation  of  this 
unexpected  discrepancy  is  definitely  warranted. 
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The  differences  between  the  generalized  Sommerfeld 
integrals  we  have  put  forth  and  those  given  by  C.T.  Tai  have 
been  discussed.  It  is  interesting  to  note  that  should  we 
use  C.T.  Tai's  formulas  in  the  unimoment  method,  the 
numerical  objective  would  be  to  find  f(X)  in  the  integrals, 
such  as. 
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2 

2 


-u 

+u 


1 
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Jm(Ap) 


-u, (2d-z) 
e  ^  dX 


(52) 


And,  using  the  formulas  we  derived,  the  numerical  objective 
is  to  find  An  in  the' summation  of  the  type 
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x  pm 

k  u2  n 


U_2  )  V~U1 

K  k2J  u2+Ul 


Jm(XC) 


-u_ (2d-z) 

2  dX 


(53) 


which  is  numerically  more  preferable  than  (52) . 


All  the  formulas  presented  have  been  numerically 
tested  and  found  to  be  applicable  to  solve  scattering  by 
buried  obstacles. 
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Figure  la.  The  terminating  surface  of  the  finite 

element  method  using  Tai's  modal  fields, 
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Figure  lb.  The  desired  terminating  surface  of  the 

finite  element  method  using  proper  modal 
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Figure  2.  Branch  Cuts  and  Permissible  Paths  of 
Integration. 
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Figure  3.  Coordinates  in  the  Meridional  Plane 
($=0  and  $=T0  • 
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H-Y  incidence 
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(a)  H-Y  incidence 


The  square  on  the  ground  plane  in  which  the  3D 
and  contour  plots  of  the  scattered  electric 
fields  will  he  shown 


Figure  8.  3-D  and  contour  plots  os  scattered  E-field 

amplitude  on  the  earth  surface 
(0.  =  45°,  frequency  =  700  MHz) 
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APPENDIX  A.  Proof  of  Recurrence  Relation  1 


a  m  m  dhrl2^  (kr*  m 

3^[hm  (kr)  VCOS  9)3  =  “dkl -  Pm(cOS  9)  Sin  6 


1  .  (2)  ,,  .  n  dPm(cos9) 

-j —  h'  '  (kr)  sm  6  cos  9  m 
kr  m 


"  Ci&hm2){kr)  ‘  Okr>3 


d  cose 

1  pm+l 
(2m+l)  m+1 


(cos  0) 


1  ,  (2)  ,,  .  ,  m  „m+ 1  .  . .  in  _m .  - .  . 

-i —  h  kr)  •  [■=— rr  P  ,  ,  (cos  9)  -  — r — —  P  (cos  0)  ] 

kr  m  2m+l  m+1  sm  0  m 

-  nn  hm2’  (kr>  Pm(cos  6>  -  72m+lT  hm^(kr>  Ocos  9> 


Recurrence  Relation  1  is  obtained  by  multiplying  both  sides  of 
the  above  equation  by  -(2m+l). 


APPENDIX  B.  Proof  of  Recurrence  Relation  2 


a  o\  m  dh<23  (kr) 

3kithn  (kr)  Pn(c0S  0)3  = - Oki -  cos  6  Vcos  9> 
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Kr  n  n  n+i 


=  -h^23  (kr)  •  cos  9  Pm(cos  9)  +  - P+-m-3-  h^23  (kr)  pm  ,  (cos  0) 
n+i  n  xr  n  n-i 
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